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Abstract 

The Klein-Gordon equation in the presence of a strong electric field, taking the form of the Mathieu equation, is studied. A novel 
analytical solution is derived for particles whose asymptotic energy is much lower or much higher than the electromagnetic field 
amplitude. The condition for which the new solution recovers the familiar Volkov wavefunction naturally follows. When not 
satisfied, significant deviation from the Volkov wavefunction is demonstrated. The new condition is shown to differ by orders 
of magnitudes from the commonly used one. As this equation describes (neglecting spin effects) the emission processes and the 
particle motion in Quantum Electrodynamics (QED) cascades, our results suggest that the standard theoretical approach towards 
this phenomenon should be revised. 


1. Introduction 


At present days, several laser infrastructures with ex¬ 
pected intensity of 10 24 - 10 25 W/cm 2 are under construction 
worldwide!]] [2 , 3; |4| . The experimental availability of such 
intense field sources creates exciting opportunities in many re¬ 
search fields 0, such as QED in the presence of strong fields 
|6l , Schwinger mechanism (7] El 0 [TUI [Tl], Unruh radiation 
H2 na, novel fast ignition schemes Humana, particles 
acceleration IfTTl . high harmonics generation lfT8l as well as nu¬ 
clear physics and the search for dark matter candidates mm. 

The fundamental physics underlying all these scientific ap¬ 
plications is the interaction of an intense electromagnetic (from 
now on we shall use the initials EM) field with an electron. The 
nature of the interaction is determined by the normalized field 
amplitude f = ea/m and the quantum parameter 


V = 4 (!) 

nr ' 

The electron mass and charge are denoted by m, e respectively, 
and its asymptotic momentum (i.e. the momentum in the ab¬ 
sence of the EM field) is p M = (po,Pi, Pi, Ps)- The amplitude 
of the vector potential is a — V-A 2 and the EM field tensor 
is F^ v = d^A v -d v A M . Natural units are used, namely h = c = 1, 
where h is the reduced Planck constant and c is the speed of 
light. If both £ and x are larger than one, the electron dynamics 
is both quantum and non-linear. The appropriate framework is 
the strong-Field QED ll2TI . Its basic principle is the inclusion of 
the term corresponding to the interaction with the classical laser 
field into the free part of the Lagrangian. As a consequence, the 
unpertubed states appearing in the cross section calculation are 
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no longer free waves. Instead, they are represented by the solu¬ 
tions of the quantum equation of motion in the presence of the 
EM field. These wavefunctions were obtained by Volkov EL 
provided that the EM field propagates in vacuum. 

Employing this approach, the properties of QED in the non- 
perturbative regime were thoroughly investigated through the 

years E3 El E3E3E21 EMESES ED El IS]. The lowest 

order strong-field processes are the non linear Compton scat¬ 
tering, where an electron interacts with the laser photons and 
emits a hard photon and the non-linear Breit-Wheeler process, 
where a photon decayes into a positron-electron pair in the pres¬ 
ence of the EM field ll23l . A sequential series of these pro¬ 
cesses, called ”QED cascade” is followed by a rapid formation 
of a QED plasma whose ingredients are electrons, positrons and 
gamma photons El. Besides its fundamental significance, this 
phenomena attracts scientific attention as a possible laboratory 
astrophysics settings E2 and as a potential gamma ray source 
j36l . Furthermore, it was suggested that spontaneous cascades 
impose a limit (of about 10 25 W/cm 2 ) on the achievable laser 
intensity (37]]. 

2. The physical scenario 

The most favorable EM field configuration for the generation 
of a QED cascade is a rotating electric field (38l . It can be re¬ 
alized experimentally in the anti-nodes of counter propagating 
laser beams. The kinetic calculation of the cascade formation 
is based on a Monte Carlo technique describing the quantum 
processes mentioned above, integrated with a Particle-In-Cell 
code taking into account the collective EM field influence on 
the classical motion of the electrons lf38l [39l l40l f4ll l42l l43l l44l . 
The quantum rates were obtained with the Volkov wavefunc¬ 
tions If23l . though not formally adequate for a rotating electric 
field. 
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The underlying assumption (see p. 3 in Il23l and p. 454 
in Ii45l ) justifying this technique is that as long as the nor¬ 
malized EM field invariants, T = e 2 F pv F llv /(4;w 4 ) and Q = 
e 2 e a j 3 pv F a ^F pv /(4m 4 ), are negligible with respect to 1 and x 2 , 
namely 

T,g« ir ( 2 ) 

the particle dynamic is well described by the Volkov wavefunc- 
tion. The symbol e a p pv above stands for the Levi-Civita ten¬ 
sor. The aim of this publication is twofold. First, we explicitly 
demonstrate that this long-believed assumption is incorrect. For 
this purpose we derive a novel solution to the quantum equation 
of motion in the presence of a rotating electric field. The com¬ 
parison between the Volkov wavefunction and our new solution 
in the relevant parameters range reveals an overwhelming devi¬ 
ation. Consequently, the common theoretical treatment towards 
the QED cascades has to be revised. Second, the new solution 
derived below suggests a natural way to obtain the emission 
rates in the presence of a rotating electric field for the sake of 
QED cascades calculations. 

3. The Governing Equation 

The particles of interest in the context of QED cascades are 
fermions (electrons and positrons). Therefore, their wavefunc¬ 
tion obeys the Dirac equation. However, the spin effect may be 
neglected in the typical cascades conditions ll46l . Hence, for 
the sake of simplicity we shall treat them as scalars. Thus, the 
free equation of motion of the particle wavefunction (D is the 
familiar Klein-Gordon equation in the presence of a EM field. 

[-<9 2 - 2 ie(A ■ d ) + e 2 A 1 - m 2 ] 0> = 0. (3) 

The center dot stands for Lorentz contraction. The EM field 
depends upon the spatial and temporal coordinates through <p = 
k ■ x where k p = ( cop , 0,0, k z ) is the wave vector. We assume a 
circularly polarized field 


of ultra-relativisitc intensity £ » 1, the quantum emission pro¬ 
cesses take place on a length scale much shorter than the laser 
wavelength If23l . meaning that the spatial dependence may in¬ 
deed be ignored. Notice that this dispersion also describes an 
EM wave propagating through plasma, with m p h analogous to 
the plasma frequency of a classical plasma wave Ez). A rigor¬ 
ous Lagrangian formulation of this plasma wave may be found 
elsewhere 11481 . 

The wavefucntion is characterized by the asympototic mo¬ 
mentum p, namely <f> ; , = e~ ,px F(<p). The substitution of this 
ansatz into ([3]) yields 


m p h 2 F" - 2ie(k ■ p)F'+ 

[i e 1 a 2 + 2 ea(p ■ e)e'^ + 2ea(p ■ F — 0. (5) 


Equation 0 was considered by several authors. In the follow¬ 
ing the main approaches are briefly reviewed. If (p ■ e) vanishes 
the quantum equation of motion may be solved analytically 1491 
even with the spin term included (i.e. Dirac equation). A so¬ 
lution for a discrete set of momentum values was derived in 
123 ED- If the EM wave dispersion relation is assumed to be 
vacuum-like, i.e. m p h vanishes, then 0 reduces to a 1 st order 
equation and therefore admits an exact solution - the familiar 
Volkov wavefunction |[22l [451 . As a result, a perturbative ex¬ 
pansion with resepect to m p h may yield an approximate solu¬ 
tion. Since the small parameter multiplies the highest deriva¬ 
tive, singular perturbation techniques are required l52l . The 
most appropriate for our problem is the familiar WKB, which 
was applied by Gam. 

A different approach was taken by II55U56I . showing that ([5]> 
is equivalent to the familiar Mathieu equation. With the aid 


of the transformation F = y(z) exp 


• k‘P m 

i~r<p 


, one can prove the 


_ j _ 

equivalence of (pb to the Mathieu equation [551 


y" + [A-2q cos 2z] y - 0 (6) 


A(<P) = a(<f>) ( ee ilp + (4) 

where the polarization vectors are e = (<?i - zea)/ V2 and * de¬ 
notes complex conjugate. a((f>) is a slowly-varying amplitude 
vanishing at <p —» +oo. However, in the following this enve¬ 
lope is assumed to be slow enough so that a 2 is approximately 
constant. 

The most general dispersion relation of the EM field is 

_>2 

massive-like k 2 — ur L - k = m 2 ph , where m p h is the effective 
mass of the EM wave photons. In the wave frame of reference, 
k z = 0 and the particle experiences a rotating electric field with 
frequency zzz p /,. On the other hand, it is well known that a stand¬ 
ing wave formed by two counter-propagating beams takes the 
form A oc cos (k z z) cos (tot). Consequently, a particle located in 
the vicinity of the anti-nodes, where cos k z z~ 1 + 0(k 2 z 2 ), also 
experiences a rotating electric field. As a result, Eqs. 00 
enable us to study the dynamics of an electron placed in the 
anti-nodes of the standing laser wave, where the wave frame of 
the former coincides with the lab frame of the latter. It is justi¬ 
fied to restrict our discussion to the anti-nodes since in the case 


where the following relations are used 


A = 


4 

m 2 . 

ph 


(k • PF , ,7 

- -z -h yea) 

m , 

ph 


q = 


-ea\p ■ e\ (7) 


ph 


and 

z = [cp + M/2, p-e=\p-e\e i *°. (8) 

The tag symbol stands for derivative with respect to z. 

In this paper, we construct an approximate solution to 01 
employing a novel mathematical technique. We shall prove that 
under the condition q/A <sc 1 it is equivalent to an effective 1" 
order equation and thus is easy to solve. 


4. The Novel Solution 

We start by reviewing the standard way to solve numerically 
the Mathieu equation ll57l 58 j. This formalism will prove use¬ 
ful for the sake of our novel derivation. Eq. 0 bears an ap¬ 
parent similiarity with the equation of motion of an electron in 
a crystal, if A p {cf>) is replaced by the periodic potential. Hence, 


2 








the Floquet theory (analogous to Bloch theorem in solid state) 
is applicable. Accordingly, the solution may be expressed as 
y = P(z)e ll ‘ 7 where P(z) is a periodic function and is called 
the characteristic exponential. Due to its periodicity P{z) may 
be expanded in a Fourier series. 


y{z) = e^Y J C 2 ne 1,zn (9) 

n 


where c 2 „ are the spectral coefficients. Substituting ([9]i into ([ 6 | 
yields 

VlnCln — ('2n+2 + Gin—2 ( 10 ) 


where the notation of 


V 2 n = 


is adopted, i.e. 
A - ( 2 n + /j ) 1 

q 


( 11 ) 


Dividing (|T()|) by c 2n - 2 we have 


Now let us examine the following 1 st order equation. 


y+i 


Vd —cos ( 2 z) 

Vd 


y = o. 


( 18 ) 


One can verify that substituting the Floquet ansatz (j9| into (18 1 
yields the recursion relation ( 101 with ji — - Vd and V 2n identi¬ 
cal to (17 1 . It implies that as long as (16 1 holds, the original 2 nd 
order equation is equivalent to an effective I s ' order one. The 
effective equation (p~8]> admits the straightforward solution 


y = exp 


-i yfAz 


+ l ~ 


2 Vd 


sin( 2 g) 


(19) 


The analytical form of © provides us with an expression 
for the coefficients. Employing the identity | 


e 


iX sin (p _ 


= 2 J n (X)e^ 


( 20 ) 


VinG 2n - G 2n G 2n+2 + 1 


(12) where J n is the Bessel function, we deduce that 


where G 2n = c 2 „/c 2 „_ 2 . Hence, G 2n is given by 


G 2 „ - 


1 

V 2n - G 2n+2 


(13) 


It can be shown that for sufficiently large n, the dominant term 
in the denominator of © is V 2n . Since V 2 „ is growing poly- 
nomially in n, it follows that c 2n decays rapidly in n. Conse¬ 
quently, the series ([9} may be trancated at a certain index de¬ 
noted by n*. Hence, G 2 „> = 1 / V 2n > , and the lower harmonics 
are obtained through the recursion relation ( fL3] >. Similiarly, for 
n < 0 we have 


H - 2 „ - 


1 

V- 2n - 2 - //_ 2 „_ 2 


(14) 


Cln = Jn (X) , X=-^~. (21) 

2 Vd 

The Bessel function vanishes for X - 0 and reaches its first 
maximum, for n » 1, at X ~ n + 0(n ^ 3 ) Il58i . The rise to 
the peak is extremely rapid, so that J n (X) is practically zero for 
X < n. In order to establish this statement, we take advatage of 
the expansion of J n (X ) for X < n (p. 250 in {591 , leading term 
only). 


J n (X) » 


tanh a(X) 1 1 , \ 

-A |/3 —n tanh a(X) exp[/j(tanh a(X)+ 

n^/3 \3 / 


1 

3 


tanh 3 a(X) - ar(^ - ))] 


( 22 ) 


where // _ 2 „ = c_ 2 „_ 2 /c^ 2 „. The solution consistency requires 
the ratios Hq, G 0 to be related by 

HqGq - 1. (15) 


Iterating this condition, the characteristic exponential /j is nu¬ 
merically obtained. 

This is the starting point of our derivation. At the moment, 
the wavefunction spectral width n* is unknown. Let us assume 
that it obeys 

2 if « /j. (16) 


Therefore, the quadratic term in n 1 appearing in (11 1 is negligi¬ 
ble and V 2n is anti-symmetric with respect to n. As a result, one 
can verify that /j — ± Vd yields Vq = 0 and satisfies (15 i. In the 
following we shall obtain an explicit formula for n* and thus 
determine the validity range of our approximation. Substituting 
Iu - - Vd in the expression for V 2n leads to 


V 2n 


4 n Vd 


(17) 


It should be mentioned that we do not consider the solution that 
corresponds to the positive root v = Vd for physical reasons, as 
shall be discussed below. 


where K\p is the modified Bessel function and a new quantity 
is introduced 



(23) 


We are interested with X « n, namely a <s 1. Hence, we Taylor 
expand the exponent argument 


1 Ol 

tanh a + - tanh 3 a - a « — — 


(24) 


and employ the asymptotic expression of K \/3 for a small argu¬ 
ment u 


K\/ 2 (u) ~ C\u 1 ^ 


(25) 


where C\ is an insignificant constant. Substituting ( |25| l and (24 1 
into ( 221 we get 


J n {X) « C 2 n 1/3 exp 



(26) 


where C 2 is another constant. Let us write J„ in terms of AX = 
n — X instead of a. 


a « tanh a * 




(27) 
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Finally, we have 


J n (n - AX) « C 2 «' 1/3 


exp 


77 

(2AX\ 5/2 ~ 

5 

l « ) 


(28) 


That is to say, the decay is exponential and the width is negligi¬ 
ble, namely AX/n oc n~ 2 ' 5 . As a result, ploting the coefficients 
c 2n (X) as a function of n for a given X one observes a rapid 
decay for n > X, as demonstrated in the results section below. 
Consequently, the spectral width of the solution © is 


2 n 


\q\ 

yfX 


(29) 


Comparing it with (16 1 we obtain the validity condition, which 
takes the following form 


6 = 2 n/p 


«c 1. 


(30) 


The exact solution of ([6]) is known to exhibit instabilities in 
certain domains of the (q. A) space. Moreover, ( fTT) implies that 
the spectral distribution is essentialy non symmetric due to the 
term proportional to n 2 . However, the expression ( fl9] l, being 
a solution of a 1 st order equation, is spectrally symmetric and 
always stable. Subsequently, we deduce that as long as © 
holds the solution lies in the stable region and the spectral non- 
symmetry is negligible. 

5. Physical Interpretation 

Let us write the final solution to |3]i with the aid of The 
wavefunction <1^, is related to y(z) by 


® P (x) = y(z) exp 


-ip ■ x + 


Ak'P) 


ph 


We define 


_ p (k-p) 

v = 9 + — r~ 
z m . 

ph 


(31) 


(32) 


and use the expressions (|7][8]i relating q, A and z to the physical 
quantities of the problem. 

4> p (x) = exp \-i(p + vk) ■ x - in sin (k ■ x + 0 O )]. (33) 

The phase cpo is debited in ([8]) and n*,v read 


eci\p ■ e\ k-p 

y= -w-(D- 1) 


(k ■ p)Sl’ 

where a new quantity is introduced 


(34) 


t 7 

ph 




(35) 


In order to clarify the physical meaning of n*, v we eploit again 
the Fluquet representation of the wavefucntion 

<X> p (x) = e ~ i(p+vk) ' x ^ c 2n e in(kx) (36) 


where the coefficients, in accordance with ([20]), take the form 

c 2n = J n (n) ■ (37) 

The wavefunction ( |36| is a superposition of free waves with 
momenta p p + (v + n)k p weighted according to c 2 „. Each wave 
may be regarded as an electron carrying v+n laser photons. The 
quantity as was shown in the previous section, corresponds 
to the width of the spectral distribution, and v may be regarded 
as its center. In other words, it is the average number of laser 
photons carried by the electron, leading to the debnition of the 
quasi momentum 

Q m = Pn + kpV. (38) 

The effective mass is associated with the quasi-momentum 
through 




i + (—)". 
\mi 


(39) 


Interestingly, the effective mass is the same as in the Volkov 
case, though v is different (as discussed below). O determines 
the deviation from Volkov. This statement becomes apparent if 
we recall B3I that that the Volkov wavefunction takes the same 
form as (|33| but with 


ea\p ■ e\ 

«„ = —-—, v v = 


e 2 ci 2 
2 (k-p) 


(40) 


(k-p) ’ 
leading to the coefficients 

4 = ■/„(%)■ (41) 

Dividing ( [34] ) by ([40]) the following relations are obtained 

2 


o"' 


Q+ 1 


Vy. 


(42) 


Taking the limit m p i, —> 0 yields Q —> 1 so that the Volkov 
wavefunction is recovered. Now the neglection of the sec¬ 
ond solution corresponding to p = Vd may be comprehended. 
Should we take the positive root, the brackets in ( [34] ) would 
become (Q + 1). As a result, the limit m p h —> 0 leads to un¬ 
physical divergence. In the opposite limit, O » 1, we have 
Q as eam p h/(k ■ p) and therefore 


ea 

l^lph 


I p-e\ 

IMph 


(43) 


Notice that in this case, as opposed to the Volkov case, n* is 
independent of the EM held amplitude and v is independent of 
the particle momentum. 

We have mentioned in the introduction that according to the 
common knowledge the condition ([2] allows one to approxi¬ 
mate the wavefunction by the Volkov solution even for a non¬ 
vanishing m p i Let us write down x- explicitly. For this pur¬ 
pose the vector potential (|4]i is substituted into |I][2]i. 


T = 


I eam p h \ 2 
\ m 2 I 


X 


(ea) 2 (k ■ p) 2 + (eam ph ) 2 \p ■ e\‘ 


(44) 


Since the second term in the expression for^f 2 is allways smaller 
than the brst, and as we are interested in order of magnitudes 


(36) 


4 
























only, we have x 2 
become 


(t °y . Therefore, the conditions (j^Ji 



However, Eq. (42 1 explicitly proves that m p h may be neglected 
only if Q w 1, or equivalently 


/ eaniph \ 2 

l k-p ) 


« 1. 


(46) 


Let us plug in typical numbers for a standing wave created by 
optical laser beams with intensity I 10 2A W/cnr (as expected 
in ELI Q), namely ea/m = 10 3 as well as m p h/m = 10 6 . It 
is favorable, for this purpose, to evaluate ( [45] [46] in the wave 
framework, where k-p - pom p h . We obtain 


the WKB approximation, corresponds to the expansion of Q in 
powers of m p h- Notice that in the above derivation, as opposed 
to the WKB approach, m p h is not assumed to be small and may 
get any value as long as 6 <K 1. 

C. If the particle energy is of the same order of magnitude as 
the field amplitude ( p 0 ~ ea) there are two possibilities, de¬ 
pending on the value of p ± . On the one hand, p ± <s p a leads 
to <5 <sc 1 and Q ~ y/2, meaning that our solution is valid. On 
the other hand, p ± ~ p () , we have 6 « 1. Consequently, the 


effective equation (181 does not represent the original Mathieu 
equation, giving rise to second order behavior such as bands 
structure formation. This regime is extremely interesting and 
will be addressed in a seperate publication. 


6. Numerical Results 


r = 
r _ 

X 2 


I eanipf, \ 
V m 2 ) 


and 


Q 2 - 1 = 


/ mm ph \“ 

U • p) 

/ \ 2 
I eam p h \ 


= 10" 


Po 


= 10 b — 


(47) 


(48) 


(49) 


l k-p J \p 0 

One can see that the difference between ( |47|48] and ([49] is of 
several order of magnitudes. As a result, according to (48 49 1 


the Volkov approximation breaks only if po ~ m while our new 
condition ( |49| states that the Volkov solution is not valid in the 
range m < po < ea, which is extremely relevant for cascades 
formation. 

Let us examine 6, Q as a function of the asymptotic momen¬ 
tum for a given m p h. In terms of the physical quantites, 6 (de¬ 
fined in (|30|) takes the form 


d = 


lecinrjp ■ 


(k ■ pY + (eam ph ) 


2 ' 


(50) 


In the wave framework. 


\e- P I 


P\ +P\ = Pa_ 

2 V2' 


Finally, we have 


V2 eap±_ 
p 2 + (ea) 2 ’ 



(51) 


(52) 


There are several regimes, according to the values of po,ea and 
P±- 

A. If po <s ea we obtain 5 <K 1 and Q » 1. The meaning is 
that our solution is valid and its deviation from Volkov is sig¬ 
nificant. For the laser parameters mentioned earlier, ea/m = 
10 3 , m p h/m = 10 ( \ the value of Q lies in the range 1 - 10 3 . 

B. In the opposite case p o » ea, it follows that ^ « 1 and 
fl » 1. Namely, our solution is valid and recovers the Volkov 
wavefunction. The first correction to the Volkov solution, i.e. 


In the following, our novel analytical solution is compared 
with the Volkov wavefunction and with the exact solution. The 
goal is twofold - illustrate its deviation from Volkov as well as 
investigate its accuracy for varying values of the small parame¬ 
ter 6. As described in the previous sections, all three solutions 


may be cast in the general form (36 1 and are therefore char¬ 
acterized by v and C 2 „. These quantities are given by ([34 37 1 
for our new solution and by (40 41 1 For Volkov. As to the 
exact wavefunction case, they are obtained numerically accord¬ 


ing to (11 13 - 15 32]. Notice that adding to v any integer j 


leaves the solution ( [36] unchanged besides a shift in the distri¬ 
bution C 2 n C 2 n+ 2 j- It allows the numerical algorithm, when 
searching for v e , to be restricted to v e = v e where v e lies the 
range 0 < v e < 1. Afterwards, it may be shifted by an integer, 
v e = v e +j, so as to make the distribution centered around n = 0. 
The normalization of the various solutions was determined ac¬ 
cording to the condition 

Jd 3 x [<Vo®; - = 1 ■ (53) 


For the sake of demonstration only, the plots appearing below 
were calculated with different laser parameters than these of 
ELI mentioned in the previous section (ea/m = 10 3 , m p h/m = 
10 fi ). The reason lies in fact that these parameters result in an 
enormous number of harmonics (as one obtains by substituting 
them into ( |34] ) 

n = 10 9 —, 1 < £1 < 1000 (54) 

Q.po 

making the wavefunction graphically difficult for inspection. 
As a result, the laser parameters were chosen to be £ = 
20, m p h = m/ 10, yielding the range 


n = 2000—, 1 < O < 20. 

Qpo 


(55) 


Nevertheless, the physical effect we wish to demonstrate re¬ 
mains the same. 

We start with calculation parameters corresponding to 6 1, 


i.e. small transverse momentum p ± <s po, ea (see Eq. (52 1 ). 
We have taken p = (m/5,m/5,0), corresponding to Q ~ 19.3 


5 


































and <5 = 0.02. The spectral shape of the Volkov and the new 
wavefunctions is shown in Fig. 1. The exact wavefunction 
was calculated as well, but was not plotted in the figure as its 
deviation from our new analytical solution are extremely neg¬ 
ligible. The deviation from the Volkov wavefunction is over¬ 
whelming - the width of the new solution is ~ 20 times smaller 
than that of Volkov, in agreement with the value of Q. It should 
be stressed that each distribution is centered around a different 
v. For Volkov we have vy = 1924.5, while our solution corre¬ 
sponds to v = 189.9. 

Fig. 2 depicts the exact, analytical and Volkov wavefunc¬ 
tions for p = (5m, 5m, 0), corresponding to 6 = 0.44,Q = 2.97. 
Several interesting points stem from the comparison. First, the 
difference between the spectral width of the Volkov and the an¬ 
alytical solutions decreases with respect to the previous case. 
Second, the analytical solution is fairly close to the exact one, 
even though 6 is not negligible. Third, the spectral shape of 
the exact solution is deformed, implying that the symmetry be¬ 
tween photon emission and absorption no longer exists. The 
distributions are centered arround vy = 280.06, v = 140.3 for 
the Volkov wavefunction and the new solution respectively. For 
the exact wavefunction, however, the non-symmetric distribu¬ 
tion implies that the definition of v as the center of the distri¬ 
bution loose its meaning. Hence, we have chosen v e to be as 
close as possible to the analytical solution, in order to ease the 
comparison. The numerical calculation yielded v e = 0.97 and 
was shifted, according to the above argument, to v e = 140.97. 


7. Conclusion 


The equation of motion of a particle in a rotating electric 
field, taking the form of the Mathieu equation, was analyzed. A 
novel approximated solution was found, adequate if the asymp¬ 
totic energy of the particle (in the wave frame) is much smaller 
or much higher than the held amplitude. For the first case. 


po <s ea, the spectral width of the analytical solution (33 i was 
shown to be il times smaller than the width of a Volkov wave- 
fucntion with the same p p . For the second case, po » ea, cor¬ 
responding to an energetic electron beam colliding with a laser, 
our solution recovers the Volkov wavefunction. The differences 
between Volkov, our new solution and the exact solution in the 
intermediate regime po ~ ea were explored numerically. 

As described in the introduction, the emission rates embed¬ 
ded in standard QED cascade calculations rely upon the Volkov 
wavefunction and depend on the variables £,x- The underly¬ 
ing assumption is that the Volkov solution is applicable as long 
as 0 holds. However, according to the above analysis, devia¬ 
tions from the Volkov solution occur unless po » ea, in con¬ 
tradiction with 0 - Consequently, we argue that the emission 
processes are no longer described by the well known expres¬ 
sions obtained with the Volkov wavefunction. Moreover, they 
depend upon another parameter, taking into account the value 
of m p h- The modifed rates may be obtained with our novel so¬ 
lution presented above. The apparent similarity between the 
mathematical structure of our solution and the Volkov wave- 
function implies that the mathematical techniques exploited to 
derive the Volkov rates may be of use in our case as well. The 


predicted cross section of such a calculation is supposed to be 
of quantum nature and to deviate from Volkov provided that 
X ~ 1 and po < ea. For optical lasers (m p h/m ~ 10 6 ) it follows 
from (44 1 that ea/m ~ 10 3 (corresponding to I » 10 24 W/cm 2 ) 
is required. Consequently, our solution may be put to test with 
the next generation laser systems. 
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Figure 1: (color online). The wavefucntion spectral coefficients of the Volkov 
(blue curve) and the analytical (red curve) solutions respectively, for f = 
20, m p h = m/10, p = (0.2m, 0.2m, 0), corresponding to 8 = 0.02. 
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Figure 2: (color online). The wavefucntion spectral coefficients of the Volkov 
(blue curve), the analytical (red curve) and the exact (black curve) solutions 
respectively, for f = 20, m^ = m/10, p = (5m, 5m, 0), corresponding to 8 = 
0.44. 
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